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Study of the Finite Density State based on SU(2) Lattice QCD 
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We report our recent numerical studies on two-color QCD with Wilson fermions. First the phase structure of the 
system in (k, /i)-plane is analyzed by measuring Polyakov line, baryon number density, gluon energy density and 
Polyakov line susceptibility. Then meson and baryon (di-quark) propagators are investigated. We find that the 
vector meson mass decreases near Hc. This phenomena can be related to the low-mass lepton-pair enhancement 
observed in CERES experiment. Since our lattice is small, in order to confirm the results, we calculate both 
periodic and anti-periodic boundary condition cases. Preliminary calculations of color averaged, symmetric and 
anti-symmetric forces of qq system obtained from Polyakov line correlations, and gluon propagators are reported. 



1. Introduction 

Finite density QCD has attracted much atten- 
tion recently: we are expecting that QCD has a 
rich structure when we add a density or chemi- 
cal potential axis in parameter space. The finite 
temperature deconfinement transition line starts 
at /i = 0, which might be a cross-over accompa- 
nied by a tri-critical point at finite /z. At high 
fj, and relatively low temperature, there would be 
a new form of matter such as color superconduc- 
tive phase, or color-flavor locking phase, which 
are speculated by mainly perturbative calcula- 
tions If. 

Non-perturbative study of the finite density 
QCD by lattice simulations is highly required, but 
because of the famous difficulty, the direct Monte 
Carlo simulation of SU(3) system is very difficult. 
Many approaches have been reported, j^-^, and 
yet, no reliable algorithm is known to investigate 
color SU(3) QCD at low temperature with large 
chemical potential. Color-SU(2) is currently the 
unique model by which we can study these re- 
gions non-perturbatively. After the first trial of 
the calculation 01, there had been no work during 
many years. Because of the recent rapid progress 
of theoretical understanding of the system |p|,p |, 
however, the calculation has been activated TlC - 



In this report, we study the finite density QCD 
by Wilson fermions. We investigate the phase 
structure in (k, /i)-plane for Iwasaki improved- 
gauge action. In order to take account of the 
fermion determinant, we use an algorithm in 
which the ratio of the determinant is evaluated 
explicitly at each Metropolis update process [p^ . 

2. Phase Structure 

In Fig.|l|, we show the average of Polyakov line 
in (K,/Lt)-plane at /3 = 0.7. We evaluate also 
baryon number density, gluon energy density and 
Polyakov line susceptibility, and especially from 
the peak of the Polyakov line susceptibility we can 
determine the confinement/deconfinement phase 
transition regions. In both bounday conditions^, 
the behaviors of these quantities are essentially 
the same |l3). 

We study several cases by varying lattice size 
(4^, 4^^ X 8, 4'^ X 12), the boundary conditions and 
the number of flavor {Nj = 2,3). 

Calculations become unstable and break down 
at large chemical potential regions [Q. In such 
cases, negative determinant appears for Nf — 3. 
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^Periodic (anti-periodic) refers the periodic (anti-periodic) 
boundary condition for quark fields in spatial directions. 



Figure 1. Behavior of Polyakov line in {k, fi)- 
space. 



3. Meson and diquark propagators 

We calculate meson and di-quark propagators 
in pseudo-scalar and vector channels. We find 
that p meson mass drops as fi reaches to the 
critical region. In order to confirm this unex- 
pected result, we calculate both periodic and 
anti-periodic boundary conditions and several k's 
(k = 0.150,0.160,0.175). In Fig.|we show tt and 
p masses as a function of fi for k = 0.160 (peri- 
odic boundary condition) . Although our statistic 
is not enough to extrapolate them to the chiral 
limit, the signal of the anomalous behavior of the 
vector channel is clearly seen. 

At p = 0, because of the QCD inequality ||], 
the pseudo-scalar channel has the lightest mass, 
but there is no such a condition for the finite 
density state. Indeed there are several conjec- 
tures in literatures iQ. If our result shown in 
Fig.|| is not a special feature of the color SU(2) 
model, this is the first lattice observation of the 
drop of the vector meson mass which may play 
an important role in the phenomenological anal- 
yses on large low-mass lepton-pair enhancement 
measured in CERES |jl|l. 

Another possibility is that the vector meson 
channel is mixed with the vector diquark one, 
'0'^757uV': at high p and their masses are dis- 
tortedn. In this case, we can study "meson" and 
"diquark" mixing in color SU(2) QCD. 

^We thank Thomas Schaefer for useful discussions on 
Vi7^i/'/V'''757MV' niixing . 



Figure 2. tt and p masses as a function of p for 
K = 0.160 with periodic boundary condition 



4. Force between qq / qq 

In the previous section, we show that the vec- 
tor meson may have a peculiar behavior at finite 
density. Above pc at low temperature, the con- 
densation of qq may occur. We need tools to un- 
derstand the vacuum structure numerically. 

From Polyakov line correlations, color average 
heavy quark potential is calculated. One can 
also obtain color singlet {Vi) and triplet poten- 
tial {V3) for qq system: x Nc = 1 (B {N^ - 1) 
p6[ |. Color symmetric (T4) and anti-symmetric 
(Va) potential for qq system can be also evalu- 
ated: N,xN, = \N,{N, -f 1) ® \N,{N, - 1). 
In color SU(2), Vi = K and = Vs even at 
/i 7^ 0. Therefore if there would exist difference 
between qq and qq states, it should originate in 
Dirac structure. 

In Fig.^ we plot preliminary data of color sin- 
glet potential. The rotational invariance is good 
by virtue of the renormalization group improved 
gauge action. The singlet (triplet) potential is at- 
tractive (repulsive) and the forces become weaker 
as p increases. We fix the gauge in Lorentz gauge 
for the color-dependent force and gluon propaga- 
tors discussed in the following. 

We study gluon propagators as a tool to inves- 
tigate the finite density QCD, 

G^,u {Po ,P^,Py,z) = {TtA^ (Po ,Px,Pv,z) 

xA,{-Po,-P^,-Py,0)). (1) 

We set the momenta as {Po,Px,Py) = 
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Figure 3. Color singlet potential Vi for /i = 
0.0,0.4,0.6 and 0.7. 



(0, 2tt/Nx, 0). The electric and magnetic parts of 
the propagator, Ge = Gtt and Gm — Gyy, are ex- 
pected to behave as, Gf.[m){z) oc exp(— i?e(„-)2;). 
When the system is screened at the finite tem- 
perature and/or finite density, the screening ef- 
fects should appear as = y/ijfTW) and 
= ^(P2 + M"^), with Me and M„ being 
electric and magnetic masses of the gluons, re- 
spectively. At T = and /i = the screening is 
infinite and the propagators vanish at long dis- 
tance 0. 

We show in Fig.^ a preliminary result for 
Gm{z). As /I increases, the screening becomes 
weaker, i.e., the system shows the deconfinement 
feature, but at fx = 0.7 the propagator at long 
range seems to drop. It is an urgent future task 
to confirm the phenomena, because it could be an 
indication of a new phase at large /x. 
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Figure 4. Transverse gluon propagator for differ- 
ent ^. Lattice size is 8 x 4 x 12 x 4. 
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